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1.0  INTRODUCTION 

Liquid  and/or  solid  particulates  are  typically  formed  in  combustion  processes  by 
condensation  of  supersaturated  vapors  followed  by  freezing  as  the  particulates  cool.  The 
small  (usually  sub-micron  in  size)  unit  particles  so  formed  may  coagulate  or  agglomerate 
into  larger  liquid  or  solid  particles.  The  nature  and  size  distribution  of  combustion-produced 
particles  is  important  from  a  system  performance  standpoint  as  well  as,  many  times,  From  an 
environmental  standpoint.  Thus,  measurement  techniques  for  particulates  are  an  important 
issue.  In  many  cases,  in  situ  optical  methods  may  be  used,  but  in  others,  the  density  of  the 
particulates  or  their  incandescence  may  preclude  direct  measurements,  and  sample 
extraction  methods  must  be  used.  Under  normal  industrial  conditions  (e.g.,  smokestacks), 
representative  particle  samples  can  be  extracted  from  the  low  velocity  flow's  using  isokinetic 
sampling  procedures.  However,  when  the  particle-laden  stream  is  flowing  at  high  velocity 
(even  supersonic),  the  flow  around  the  probe  becomes  quite  complex,  and  the  nature  of  the 
sample  collected  is  directly  related  to  the  complex  flow  pattern.  This  is  the  situation 
encountered  in  the  exhaust  plumes  of  rocket  and  turbine  engines  and  in  some  wind  tunnels. 
In  this  study,  the  problems  encountered  in  sampling  the  exhaust  flow  from  a  solid  propellant 
rocket  motor  being  tested  at  simulated  high  altitude  are  examined. 

Probe  designs  which  must  be  immersed  in  the  exhaust  plumes  of  solid  propellant  rocket 
motors  must  withstand  the  high  velocity  («  6000  ft/sec),  high  temperature  (>  3,500°K) 
environment.  The  systems  are  further  constrained  by  the  physical  arrangement  and 
dimensions  of  the  test  facilities.  Thus,  isokinetic  sampling  can  rarely  be  accomplished,  and 
the  effects  of  the  flow  on  the  sample  integrity  must  be  assessed.  In  this  case,  particle 
sampling  from  a  supersonic  stream  will  produce  a  detached  shock  wave  in  front  of  the 
particle  collection  probe.  The  shock  front  creates  a  velocity  difference  between  the  free 
stream  and  the  gas  entering  the  sampling  probe,  and  a  size-dependent  preferential 
withdrawal  of  particles  may  take  place.  In  addition,  the  large  drag  forces  which  result  from 
the  particle-shock  interaction  can  cause  the  breakup  of  droplets  and  particle  agglomerates. 

The  purpose  of  this  study  is  to  identify  the  form  of  the  scaling  law-s  and  important 
dimensionless  groups  which  are  expected  to  play  a  role  in  the  collection  efficiency  of  a 
supersonic  probe.  Important  particle  parameters  which  apply  to  the  dynamics  of  particle 
motion  behind  a  normal  shock  wave  are  identified.  These  dimensionless  groups  are  related 
to  the  stagnation  conditions  or  boundary  conditions  imposed  in  most  rocket  motor¬ 
generated  supersonic  flows.  Finally,  criteria  for  agglomeration  and  droplet  or  agglomerate 
breakup  are  established  in  terms  of  the  forces  of  adhesion  and  particle  drag. 
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2.0  SUPERSONIC  PROBE 

The  simplest  probe  geometry  is  a  straight-walled  cylindrical  tube  aligned  parallel  to  the 
free  stream  as  shown  in  Fig.  1.  More  complicated  probe  geometries,  such  as  the  type 
described  in  Ref.  1,  involve  the  use  of  a  diverging  nozzle  within  the  probe  inlet.  This 
additional  feature  is  designed  to  provide  isokinetic  sampling  conditions  up  to  a  Mach 
number  of  1.2  with  a  “shock  down,r  region  downstream  from  the  probe  inlet.  In  many 
applications,  however,  such  as  particle  sampling  in  solid  rocket  motor  plumes  with 
stagnation  temperatures  of  *4,000°K  and  pressures  in  the  range  from  500  to  1,000  psia 
(Refs.  2  and  3),  the  gas  Mach  number  at  the  rocket  nozzle  exit  plane  is  about  4  or  greater, 
and  gas  velocities  are  greater  than  3000  m/sec.  Furthermore,  the  solid  particulate  content 
may  be  as  large  as  30  percent  by  weight.  To  withstand  such  environments,  ablative  probe 
designs  are  used  in  which  the  probe  frontal  area  is  blunted,  negating  the  possibility  of 
isokinetic  sampling.  Furthermore,  the  high  Mach  numbers  make  it  difficult  to  design  an  inlet 
and  internal  geometry  within  the  probe  which  will  capture  the  particles  in  an  isokinetic 
manner. 


Detached  Shock 


Figure  1.  Schematic  of  particle  probe  in  supersonic  stream. 


2.1  CHARACTERISTICS  OF  PROBE  FLOW  FIELD 

The  flow  field  near  the  inlet  of  the  sampling  probe,  as  shown  in  Fig.  1,  is  similar  to  the 
flow  field  which  exists  in  the  vicinity  of  a  supersonic  pitot  tube  (Ref.  4).  An  essential 
difference,  however,  is  that  the  particle  probe  in  the  present  case  admits  gas  into  the  inlet. 
The  effects  of  heat  transfer  and  friction  will  normally  choke  the  flow  in  the  probe  in  such  a 
way  that  the  initial  Mach  number  at  the  probe  inlet  is  subsonic  (Ref.  5)  rising  to  sonic  at  the 
probe  outlet.  In  this  case  the  probe  is  unable  to  pass  all  of  the  free-stream  flow 
corresponding  to  its  cross-sectional  area,  and  part  of  the  flow  will  spill  over  the  inlet  lip.  The 
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fraction  of  fluid  which  spills  over  will  determine  the  position  of  the  fluid  streamlines  behind 
the  shock  wave  near  the  probe  inlet;  moreover,  the  fraction  of  fluid  lost  to  the  mainstream 
will  be  independent  of  the  probe  exhaust  pressure,  provided  that  it  is  sufficiently  low. 

The  entrained  fluid  from  the  free  stream  is  indicated  by  the  envelope  shown  in  Fig.  1. 
The  flow  is  subsonic  in  the  region  bounded  by  the  detached  shock,  the  probe  walls,  and  the 
sonic  line  which  intersects  the  probe  inlet. 

Typical  fluid  properties  across  the  shock  wave  along  the  probe  centerline  are  indicated  in 
Table  1  for  sampling  near  the  nozzle  exit  of  a  solid  rocket  plume.  Across  the  shock  wave,  in 
this  application,  the  fluid  density  remains  sufficiently  low  that  the  mean  free  path  of  the  gas 
(X)  is  about  0.5  pm.  If  the  flow'  is  choked,  the  gas  Mach  number  inside  the  probe  inlet  will  be 
less  than  the  value  of  indicated  in  Table  1.  The  fluid  Reynolds  number  (ReL)  behind  the 
shock  wave  based  on  a  probe  diameter  of  2.54  cm  is  sufficiently  large  that  the  boundary 
layers  are  very  thin  and  viscous  effects  on  the  character  of  the  flow  field  are  relatively 
unimportant. 

Table  1.  Typical  Gas  and  Particle  Properties  Across  a  Probe  Shock 
in  Solid  Rocket  Plume  Sampling 


Stagnation 
Conditions 
(Ref.  3) 

Gas  Properties 

Upstream  of  Shock 

Gas  Properties 
Downstream  of  Shock 

Particle  Parameters 
Downstream  of 
Shock  (d  =  ljim) 

T0  =  3420“ K 

M|  =  4 

M2  =  0.35 

Re2  =  5.6 

Pol  =  590  psi 

V]  =  3050  m/sec 

V2  =  402  m/sec 

Mp2  =  2.3 

P02  =  27,6  psi 

C]  =  762  m/sec 

C2  =  1 148  m/sec 

Kn  s  0.5 

7  =  1.16 

Tj  =  1500°K 

T2  =  3390°K 

St*  =  0.83 

P]  =  1.5  psi 

P2  =  26  psi 

Rep  =  6.44 

Gl  =  0.21 

/o 

II 

St*/Repz/3  =  0.24 

x  10-4  gm/cm3 

x  10  4  gm/cm3 

Rl  =  2.1  x  104 

Wb**  =  20-200 

*  Probe  diameter  =  2.S4  cm 

**2  X  103  dynes/cm2  <  P/Ap  <  2  X  10®  dynes/cm2  from  Fig,  18 

The  particle  dimensionless  parameters  for  a  unit  density  sphere  of  diameter  d  =  1pm  are 
also  shown  in  Table  1 .  These  particle  parameters,  defined  in  detail  in  the  text,  are  as  follows: 
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Reynolds  Number 
Mach  Number 
Knudsen  Number 
Stakes  Number 

Weber  Number 


dynamic  forces/viscous  forces 

particle-gas  velocity  difference/local  sonic  velocity 

gas  mean  free  path/particle  diameter 

particle  stopping  distance/characteristic  length  (probe 
diameter) 

inertial  (drag)  forces/surface  binding  forces  (adhesion) 


2.2  SHOCK  DETACHMENT 


The  shock  detachment  distance  (6)  in  Fig.  1  cannot  be  predicted  without  extensive 
numerical  calculations.  The  magnitude  of  5,  however,  is  scaled  by  the  diameter  of  the  probe 
and  is  dependent  on  the  free-stream  Mach  number  Mi.  Figure  2  illustrates  the  results  of 
experimental  measurements  for  a  flat- nosed,  solid  body  of  revolution  (Ref.  4).  The  fluid 


Mj 

Figure  2.  Shock  wave  detachment  distance  data  for  body 
of  revolution  with  flat  nose. 
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entrained  into  the  probe  may  reduce  the  ratio  somewhat,  but  a  full  numerical  solution  would 
be  necessary  to  determine  the  effect.  For  the  experimental  results  of  Fig.  2,  6/L  =  0.35  for 
Mj  >  2  with  a  maximum  of  14  percent  error.  It  is  unlikely  that  the  relative  position  of  the 
gas  streamlines  will  be  affected  strongly  for  larger  Mach  numbers  provided  the  fraction  of 
fluid  lost  to  the  mainstream  at  the  probe  inlet  remains  constant. 

2.3  PARTICLE  TRAJECTORIES 

Particles  have  significant  inertia  relative  to  the  surrounding  fluid  since  liquid  and  solid 
particle  densities  are  much  larger  than  that  of  the  carrier  gas.  Particles  within  the  entrained 
envelope  upstream  of  the  shock  wave  in  Fig.  1  will  either  follow  the  fluid  streamlines  if  they 
are  small,  or  project  across  streamlines  into  the  probe  inlet  if  they  are  large  enough.  If  the 
volume  of  gas  spilling  over  the  probe  lip  is  large  and  the  probe  fluid  intake  is  small,  some  of 
the  particles  initially  moving  near  the  inner  edge  of  the  entrained  envelope  will  subsequently 
make  contact  with  the  inner  walls  of  the  probe  lip.  The  radius  of  curvature  of  the  fluid 
streamlines  can  be  very  small  in  this  region  of  the  flow  field.  Particle  deposition  of  a  similar 
type  has  been  observed  in  both  axisymetric  and  two-dimensional  ambient  particle  samplers 
(Refs.  6  and  7).  More  significantly,  large  particles  outside  the  entrained  envelope  will  project 
across  fluid  streamlines  into  the  probe,  increasing  the  large  particle  concentration  of  the 
sampled  gas  over  that  of  the  free  stream.  (Ref.  8). 

3.0  PARTICLE  DYNAMICS 

Particles  suspended  in  the  supersonic  streams  of  rocket  plumes,  high-speed  wind  tunnels, 
and  shock  tubes  often  encounter  shock  waves.  When  a  particle  encounters  a  shock  front,  it 
projects  ahead  of  the  carrier  gas  moving  behind  the  discontinuity  because  of  ils  inertia  and 
the  sharp  decrease  in  gas  velocity.  This  phenomenon  subjects  the  particle  to  a  large  drag 
force  behind  the  shock  wave,  and  the  particle  motion  relaxes  eventually  to  that  of  the  carrier 
gas. 


3.1  NORMAL  SHOCK  RELATIONS 

For  simplicity,  a  perfect  gas  is  assumed  with  constant  specific  heats.  The  relationships 
between  gas  properties  across  a  normal  shock  wave  are  listed  below.  Reference  will  be  made 
to  these  expressions  in  the  calculations  that  follow.  Referring  to  Fig.  3, 

=  li  =  ~  0  M?  +  2 

62  v,  (7  +  i)  Mf 
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r,  ['-(V)M?]  [(t^Tt)  Mi  -  '] 

Tl  (7  +  l)2  M2 

2(7  -  1) 

where  M!  is  the  Mach  number  ahead  of  the  shock  and  7  is  the  constant  ratio  of  specific  heats 
for  the  gas. 


M]_  >  1 
M2<  1 
Vj  >V2 


Figure  3.  Particle  moving  behind  normal  shock  wave. 


3.2  PARTICLE  PARAMETERS  BEHIND  SHOCK 

Three  local  particle  parameters  can  be  defined  behind  a  shock  front.  The  three 
dimensionless  groups  are  local  values  of  the  particle  Mach,  Reynolds,  and  Knudsen 
numbers.  As  shown  below,  only  two  of  these  groups  are  independent  and  required  to 
determine  the  particle  drag  coefficient. 

3.2.1  Particle  Mach  Number 

The  particle  Mach  number  is  defined  as  the  speed  of  the  particle  relative  to  the  carrier  gas 
divided  by  the  local  speed  of  sound.  Particle  Mach  number  is  not  defined  elsewhere  in  the 
literature.  In  particle-shock  interactions,  the  particle  Mach  number  is  a  maximum  directly 
behind  the  shock  wave  and  decreases  to  zero  as  the  particle  slows  to  the  speed  of  the  ambient 
fluid.  In  general,  the  particle  Mach  number  can  be  written  as 
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Mp 


I  v  —  vf  I 

C2 


(3) 


where  v  and  vf  are  the  local  particle  and  gas  velocities.  Immediately  behind  the  shock  wave, 
returning  to  Fig.  3  and  the  expression  above,  one  obtains  a  maximum  value  for  Mp  of 


-  0  -  t)  Hr)  H)  <»> 

Since  the  local  speed  of  sound  is  c  =  ('yRT),/2  where  R  is  the  specific  gas  constant  and  M| 
=  V|/cj,  Eq.  (3a)  becomes 


m«  =  (<  -  ^7)  (t£)'/2  m:  (o 

Substituting  from  Eqs.  (1)  and  (2)  above,  one  has 

M  (HtHH-O  (5) 

For  weak  shocks  (Mi  -  1),  Eq.  (5)  indicates  that  the  particle  Mach  number  behind  the  shock 
wave  vanishes.  As  the  gas  Mach  number  increases,  however,  Mp  approaches  an  upper  limit. 
From  Eq.  (5), 


(6) 


In  the  case  of  a  rocket  exhaust  for  7  *=  1.16  (Ref.  3),  the  upper  limit  for  the  magnitude  of 
MP2  is  Mpmax  =  3.28.  For  air  with  7  =  1.4,  MPniw  =  1.89. 


3.2.2  Local  Reynolds  Number 


The  local  Reynolds  number  of  a  particle  is  defined  in  terms  of  the  particle  diameter  and 
its  velocity  relative  to  the  ambient  fluid.  As  is  the  case  for  the  particle  Mach  number,  the 
particle  Reynolds  number  is  a  maximum  behind  the  shock  front  and  decreases  to  zero  as  the 
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particle  equilibrates  with  the  gas  motion.  Thus,  behind  the  shock  front,  one  obtains  a  local 
Reynolds  number  in  the  form 


Re 


C2  1  v  -  vf  1  d 

n 


(7) 


Immediately  behind  the  shock,  however,  the  particle  Reynolds  number  becomes 


ReP2 


02  (V1  -  V^d 
V-i 


(7a) 


3.2.3  Particle  Knudsen  Number 

The  particle  Knudsen  number  is  defined  as  the  ratio  of  the  gas  mean  free  path  to  the 
particle  diameter,  or 


Kn  =  X/d 

where,  downstream  of  a  normal  shock  wave  (Ref.  9), 

,  .  n(  ™  V2 

\  2kbT2  ) 

Here,  m  is  the  molecular  mass  and  kb  is  Boltzmann’s  constant.  Taking  the  ratio  of  the  local 
values  of  the  particle  Mach  and  Reynolds  numbers  as  defined  by  Eqs.  (3a)  and  (7a),  one  has 

-_p  =  -  K  (10) 

Re  e2c2d 

where  the  speed  of  sound  is  c2  =  (7RT2}1/2.  Rearranging  Eq.  (10)  and  comparing  with  Eqs. 
(8)  and  (9),  one  obtains 

Thus,  only  two  of  the  three  dimensionless  particle  groups  are  independent  in  the  flow  field. 
It  should  be  noted  that  unlike  the  particle  Mach  and  Reynolds  numbers,  the  Knudsen 
number  remains  constant  downstream  of  the  shock  wave  provided  the  gas  properties  do  not 
change  appreciably.  The  assumption  of  constant  gas  properties  is  a  reasonable  one  in  most 
cases. 


(8) 


(9) 
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Substituting  the  maximum  particle  Mach  number  Eq.  (6)  into  Eq.  (11)  yields  the 
inequality 


Kn  Re  £ 


(12) 


In  other  words,  the  product  of  the  local  particle  Knudsen  and  Reynolds  numbers  must  be 
less  than  a  constant  which  depends  on  the  ratio  of  specific  heats  of  the  gas. 


3.3  PARTICLE  DRAG  COEFFICIENT 


Newton’s  law  for  a  spherical  panicle  moving  relative  to  a  gas  is 

mP  T  =  "T  Ap  Cd  Vf  (v*“  v?)  (13) 

dt  2 

where  v  is  the  particle  velocity,  vr  is  the  fluid  velocity,  mp  is  the  mass  of  the  particle,  q  is 
the  fluid  density,  CD  is  the  drag  coefficient,  and  Ap  (=  rd2/4)  is  the  cross-sectional  area  of 
the  particle. 

In  general,  the  drag  coefficient  (Cp)  depends  on  local  values  of  any  two  of  the  three 
groups  given  by  the  particle  Mach,  Reynolds,  and  Knudsen  numbers  as  discussed  in  Sec. 
3.2.3.  Extensive  research  has  been  conducted  on  values  of  Cp  for  particles  moving  at  both 
supersonic  and  subsonic  speeds  covering  a  wide  range  of  particle  Knudsen  numbers  in  both 
the  free  molecular  and  continuum  regimes  (Refs.  10,  11,  and  12).  A  useful  expression 
developed  by  Crowe  (Ref.  1 1)  which  correlates  drag  data  in  supersonic  flows  and  reduces  to 
Stokes  law  in  the  continuum  is  given  below. 


CD  =  (C£  -  2)  exp[-3.07  y >'2  Mp  w(Re)/Re] 

+  f(Mp)/Tl/2  Mp]  exp(  -  Re/2Mp)  +  2  (14) 


where 


fti  w(Re)  =  2.88  [l  +  tanh  (0.33  fti  Re  -  1.92)] 


and 


f(Mp)  =  [2.3  +  1.7  (Tp/Tg)l/2]  -  2.3  tanh  (0.51  to  Mp) 
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and  (Refs.  8  and  13) 


Cl  =  (1  +  0.158  Re2'3)  (15) 

Here,  Tp/Tg  is  the  ratio  of  particle  to  gas  temperature. 

Contours  of  the  drag  coefficient  for  constant  Knudsen  number,  from  Eq.  (14),  are 
plotted  in  Fig.  4.  Values  of  the  drag  coefficient  for  particles  crossing  a  shock  follow  a 
contour  of  constant  Kn  in  Fig.  4  from  an  initial  large  particle  Re  given  by  Eq.  (7)  to 
vanishingly  small  values  when  Re  —  0.  Also  shown  in  Fig.  4  is  the  inequality  of  Eq.  (12). 
Values  of  the  drag  coefficient  for  a  given  particle  must  begin  above  the  dashed  line  to  the  left 
of  the  bifurcation  point  (located  at  Re  *  28)  or  below  the  dashed  line  to  the  right  of  Re  *» 
28.  In  addition,  for  Kn  <  <  1,  Stokes  law  (given  by  Cd  =  24/Re)  and  modifications  to 
Stokes  law  in  the  inertial  regime  are  indicated. 


Several  simple  expressions  for  the  drag  coefficient  can  be  written  in  limiting  cases,  as 
shown  below. 
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Stokes: 


24 

Re 


;  Kn  <  <  1 ,  Re  <  1 


Stokes  with  slip: 


CD  =  ;  Kn  >  1,  Re  <  1 

CRe 


where  (Ref.  14) 


[ 


C  =  1  +  2Kn  I  1.257  +  0.4  expi 


<-^>] 


Non-Stokesian: 


CD  =  (1  +  0.158  Re2'3)  ;  Kn  <  <  1,  Re  >  1 

Re 


(16) 


(17) 


(18) 


(19) 


3.4  EQUATION  OF  PARTICLE  MOTION 

Considering  particle  motion  across  a  normal  shock  wave,  neglecting  the  contribution  of 
particle  acceleration  to  the  drag,  and  including  an  arbituary  fluid  velocity,  v^,  enables  Eq. 
(13)  to  be  written  in  the  form,  -Cd/24  Re(37rpt2d)(v  -  vf)  where  Re  is  the  local 
Reynolds  number  of  the  particle  and  m  is  the  fluid  viscosity  behind  the  shock.  Moreover,  we 
assume  as  before  that  only  small  changes  in  absolute  gas  temperature  occur  in  the  flow  field 
behind  the  shock  wave.  The  local  particle  Reynolds  number  appearing  in  Eq.  (20)  as  defined 
by  Eq.  (7)  is  conveniently  expressed  in  terms  of  a  reference  particle  Reynolds  number  (Rep) 
in  Eq.  (21). 


mP  ?  - rr  Cd  Re  «d(v  ”  Vf) 

at  24 


(20) 


where 


Re  =  Rep 


Rep 


62  d  V| 

n 


(21) 


(22) 
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and 

I  V  -  Vf  I  =  [(ft  -  ur)2  +  (v  -  vr)2  +  (w  -  wt-)2],/2 

for  fluid  and  particle  velocities  u,  v,  and  w  in  three  orthogonal  directions  x,  y,  and  z.  Rep  in 
Eq.  (22)  is  defined  in  terms  of  the  free-stream  gas  velocity  upstream  of  the  shock,  and  v+ 
and  vr  are  the  normalized  particle  and  fluid  velocities,  respectively  (Ref.  8). 

If  Eq.  (20)  is  divided  by  the  particle  mass  (mp),  the  initial  particle  velocity  (vj),  and  a 
characteristic  length  (L)  corresponding  to  the  probe  diameter,  and  if  changes  in  fluid 
properties  behind  the  shock  wave  are  overlooked,  then  the  nondimensional  form  of  Eq.  (20) 
becomes 


d  v 
~d¥~ 


M 5)  (±) 


subject  to  the  boundary  conditions 


(23) 


0  =  0:  |v~  |  =  1,  |  vf+|  =  -— 

vi 


(24) 


Here,  0  =  vjt/L,  CD  =  Co(Re,  Kn),  and  St  is  the  particle  Stokes  number  representing  the 
ratio  of  the  particle  stopping  distance  to  the  characteristic  length,  or 


St  = 

18  H2  L 


(25) 


It  is  also  important  to  note  in  Eq.  (23)  that  Vf  is  determined  by  the  free-stream  Mach 
number  (Mi)  and  the  fraction  of  fluid  geometrically  incident  on  the  probe  ahead  of  the 
shock  which  subsequently  spills  over  the  inlet  lip. 


Equation  (23)  has  been  written  to  include  particle  motion  in  more  than  one  dimension. 
The  particle  trajectories  determined  from  Eq.  (23),  in  any  case,  are  defined  for  a  given  flow 
field  vf+by  the  dimensionless  particle  groups  Rep,  St,  and  Kn  given  by  Eqs.  (11),  (22),  and 
(25). 


It  is  now  instructive  to  look  at  simple  approximations  to  the  dimensionless  quantity 
CoRep/St  appearing  in  Eq.  (23).  Introducing  Eqs.  (16)  through  (19)  for  CD,  one  is  able  to 
construct  Table  2,  listing  the  important  groups  which  determine  the  particle  trajectories 
from  a  given  gas  flow  field. 
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Table  2.  Dimensionless  Groups  that  Determine  Particle 
Trajectories  for  a  Given  Gas  Flow  Field 


Group 

Limitations 

St* 

Kn  <  <  1,  Rep  <  1 

C  St” 

Kn  a  1,  Rep  <  1 

St,  Rep 

Kn  <  <  1 ,  Rep  >  1 

St/Re*'3 

Kn  <  <  1,  Rep  >  >  1 

‘increase  in  Stokes  number  indicates  particles  with 
greater  inertia. 

**C  =  1  +  2Kn  [1.257  +  0.4  exp(- 0.55/ Kn)] 


4.0  SIMILITUDE 

It  is  necessary  to  understand  the  form  of  the  scaling  law  in  the  collection  of  particles  with 
a  supersonic  probe  for  the  purpose  of  laboratory  testing  and  interpreting  experimental  data. 
The  collection  efficiency  will  depend,  in  general,  on  the  properties  of  the  gas  flow  field  and 
the  trajectories  of  the  particles. 

4.1  SCALING  LAW  FOR  PROBE  COLLECTION  EFFICIENCY 

Considerable  work  has  been  conducted  on  the  problem  of  isokinetic  sampling  in 
subsonic  incompressible  flows  (Ref.  15).  Very  little  research  has  been  conducted  on  probe 
characteristics  in  compressible  flows.  Dimensional  arguments  suggest  that  the  probe 
collection  efficiency  represented  by  the  ratio  of  particle  concentration  (defined  as  particles 
per  unit  mass  of  carrier  gas  for  compressible  flows)  to  that  in  the  mainstream  must  be  a 
function  of  the  gas  sample  rate  in  the  probe  in  addition  to  groups  characterizing  the  gas  flow 
field  and  particle  motion.  Assuming  that  viscous  effects  are  negligible  in  the  gas,  one  obtains 
an  expression  for  particle  collection  efficiency  of  the  form 

St-  Kn>  Rep)  (26) 

nm  V  Qs  ' 

where  ns  is  the  number  of  particles  per  unit  mass  in  the  probe  sample  line,  nm  is  the  number 
of  particles  per  unit  mass  in  the  mainstream  upwind  of  the  detached  shock,  qm  is  the  mass 
flow  rate  of  gas  in  the  mainstream  in  the  probe  cross-sectional  area  A,  and  q,  is  the  mass 
flow  rate  of  gas  in  the  probe  sample  line. 
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For  gas  Mach  numbers  M)  >  2,  the  dimensions  of  the  flow  field  back  of  the  shock  wave 
near  the  probe  inlet  remain  relatively  constant,  as  indicated  in  Fig.  2.  If  the  particles  are 
limited  by  one  of  the  conditions  of  Table  2  (e.g.,  Stokesian  with  Rep  <  1),  the  probe 
collection  efficiency  reduces  to 

(27, 

where  M[  may  be  of  secondary  importance.  An  example  of  the  form  of  the  collection 
efficiency  is  shown  in  Fig.  5.  It  should  be  noted  in  Fig.  5  that  small  particles  (St  <  <  1) 
fallow  the  fluid  streamlines  and  n</nm  —  1.  Large  particles  (St  >  >  1),  however,  project 
straight  into  the  probe  inlet  when  conservation  of  particle  mass  suggests  nsqs  =  nmqm.  These 
limits  are  indicated  in  Fig.  5. 


Figure  5.  Example  of  probe  collection  efficiency 
for  Stokesian  particles. 


The  shape  of  the  efficiency  curve,  as  shown  in  Fig.  5,  will  be  sigmoidal,  which  is  typical 
of  the  performance  characteristics  of  aerosol  samplers  (Ref.  16).  Moreover,  St',  indicated  in 
Fig.  5,  corresponds  to  the  critical  Stokes  parameter  for  which  the  collection  efficiency  is 
intermediate  between  its  maximum  and  minimum  values.  For  particle  projection  into  a  void, 
St’  *  0(1).  Values  of  St*  are  useful  for  a  given  flow  field,  since  one  can  determine  the  critical 
particle  diameter  from  its  definition  given  by  Eq.  (25).  Particles  smaller  than  d*  follow  the 
flow  field,  while  larger  particles  project  into  the  probe  from  outside  the  entrained  envtiope. 
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In  Eqs.  (26)  and  (27),  qs  is  the  mass  flow  rate  in  the  sample  line  which  depends  on  friction 
losses  (i.e.,  friction  factor  and  probe  length),  heat  transfer  from  the  line,  and  gas  Mach 
number  (M|).  The  quantity  q,  is  difficult  to  predict  but  can  be  measured  experimentally.  The 
mass  flow  rate  of  the  free  stream  (g™),  however,  can  be  determined  from  gas  properties  in 
many  important  types  of  supersonic  flows  as  shown  in  Section  4.2.  These  considerations 
lead  directly  to  the  design  of  experiments  to  measure  the  particle  collection  efficiency  of 
different  probe  geometries. 

4.2  NOZZLE-GENERATED  SUPERSONIC  FLOWS 

The  basic  aerodynamic  element  used  to  obtain  prescribed  supersonic  flows  is  the 
converging-diverging  channel  of  the  type  shown  in  Fig.  6.  The  nozzle  is  supplied  with  gas  at 
a  high  pressure  (the  stagnation  pressure)  at  the  inlet  where  the  gas  velocity  is  small.  Along 
with  pressure,  a  stagnation  temperature  of  the  gas  is  also  prescribed.  Provided  the  exhaust 
pressure  is  sufficiently  low,  sonic  conditions  exist  in  the  nozzle  thToat,  and  the  gas  Mach 
number  at  any  position  along  the  axis  of  the  nozzle  for  a  given  ratio  of  specific  heats  is 
determined  by  the  ratio  of  the  local  cross-sectional  area  to  that  of  the  throat.  The  same  basic 
configuration  exists  in  the  nozzle  of  a  solid  fuel  rocket  motor. 


aE 


Figure  6.  Schematic  of  Laval  nozzle. 

The  nozzle  shown  in  Fig.  6  is  a  useful  tool  for  the  testing  of  probes  to  be  used  in 
supersonic  flow.  The  prescribed  boundary  conditions  in  these  tests  are  the  stagnation 
conditions  of  the  gas  and  the  nozzle  geometry.  Prediction  of  particle  parameters  can  be 
made  in  terms  of  these  boundary  conditions  as  demonstrated  in  Section  4.4.  These  results 
are  useful  in  the  scaling  of  probe  experiments  and  for  the  comparison  of  probe  performance 
with  field  tests  (e.g.,  rocket  testing). 
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If  the  nozzle  shown  in  Fig.  6  is  designed  to  function  without  significant  separation  along 
the  inside  walls,  the  flow  is  approximately  isentropic,  and  gas  properties  immediately 
upstream  of  a  shock  are  related  to  stagnation  conditions  by  the  following  expressions: 

^  =  1  +  (28) 

and 

17  =  [1  +  (-JrL)Mi]TT  <29> 


The  quantity  qm,  the  mass  flow  rate  in  the  free-stream  incident  on  the  probe  area  A,  can 
be  determined  with  Eqs.  (28)  and  (29).  Since 


qm  -  Pj  v,  A  -  (^)  (c|)  (c‘)  eoCoA 

one  has 

— S—  .  M,  (-ILY'2 

Go  A  ^RTo  V  P0  '  V  T0  / 

Substituting  from  Eqs.  28  and  29,  one  obtains 


Qm 

Co  A  yjR  T0 


Mj 

[i  +  (^t1) 


(30) 


(31) 


(32) 


This  expression  is  plotted  in  Fig.  7. 


4.3  GAS  VISCOSITY 


All  of  the  particle  parameters  discussed  in  Section  3,  with  the  exception  of  the  particle 
Mach  number  (Mp),  involve  the  gas  viscosity.  The  gas  viscosity  depends  strongly  on 
temperature  and  is  a  weak  function  of  pressure  at  very  low  temperatures.  The  Sutherland 
formula  (Ref.  9)  gives  accurate  results  for  common  gases  such  as  nitrogen  and  oxygen, 


/tV2 

[  Tr  —  Tel 

w 

L  T  +  Te  J 

(33) 
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where  T,  is  a  reference  temperature  (Tr  =  298. 2°K).  The  quantity  /ir  =  fif(Tr)  is  the  viscosity 
of  the  gas  at  the  reference  temperature  and  T©  is  Sutherland’s  constant.  Values  of  T© 
determined  empirically  for  a  number  of  common  gases  are  shown  in  Table  3,  below. 

Table  3.  Sutherland's  Constant* 


Gas  Te  (aK) 


Air 

111.3 

C02 

239.7 

CO 

100.0 

He 

80.3 

h2 

72.2 

n2 

110.6 

02 

127.0 

*From  Ref. 
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For  simplicity,  Eq.  (33)  can  be  rewritten  in  the  form 

bT^ 

H  ~  — - — 

T  +  T© 


(34) 


where 


(35) 


In  the  case  of  air,  for  example,  substituting  Tr  =  298. 2°K,  fa  =  1.8  x  10~4  gm/cm-sec  and 
T©  =  1 11. 3°K,  one  obtains  a  typical  value  for  b  =  1.43  x  10~5  gm/cm-sec-(°K)1/2.  These 
results  will  be  used  in  the  sections  below. 


4.4  PARTICLE  PARAMETERS  AND  STAGNATION  CONDITIONS 


The  particle  parameters  which  appear  in  the  equation  of  motion,  Eq.  (23),  depend  on  gas 
properties  behind  a  normal  shock  wave.  The  gas  properties  are  related  to  the  gas  stagnation 
conditions,  while  the  Mach  number  (Mi)  upstream  of  the  shock  is  related  to  the  nozzle 
design.  Thus,  expressions  can  be  developed  for  the  important  particle  parameters  in  terms  of 
the  boundary  conditions  imposed  in  the  nozzle  flow. 
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4.4.1  Stokes  Number 


The  particle  Stokes  number  is  defined  as  the  ratio  of  the  particle  stopping  distance  to  a 
characteristic  length  in  the  flow  field.  Multiplying  and  dividing  by  values  of  the  gas 
properties  across  the  shock  wave  and  in  the  stagnation  region,  one  Finds  that  Eq.  (25) 
becomes 


St 


eP  d2vi  _  ePd2 
18 n  L  _  18L 


GS)  [(£)  (Sr)  Gf)  (a] 


(36) 


Introducing  the  definition  of  viscosity  from  Eq.  (34)  and  noting  that  V]/C[  =  M]  and 
Ci/co  -  (Ti/T0)1/2,  we  can  now  write  Eq.  (36)  in  the  form, 


where 


*  -  (£)  (£) '■  (-) 


k  =  b/R1'2 


and 

•<+)(*)] 


(37) 


(38) 


It  is  interesting  to  note  that  Eq.  (37)  is  the  product  of  the  particle  Stokes  number  based  on 
the  speed  of  sound  of  the  gas  at  standard  conditions  (i.e.,  a  particle  moving  at  Mach  1  at  T 
=  298. 2°K)  and  the  dimensionless  function  fj(Mi).  Note  also  that  fi  (Mi)  *=  1  as  Mi  — ■  1. 


Numerical  solutions  of  Eq.  (38)  were  determined  by  introducing  Eqs.  (1),  (2),  and  (28). 
Contours  of  fi  (Mi)  are  shown  in  Fig.  8  for  a  ratio  of  specific  heats  y  =  1.4  for  air  and  y 
=  1.16,  typical  of  solid  rocket  exhausts.  The  results  shown  in  Fig.  8  indicate  that  for 
constant  particle  and  gas  stagnation  properties,  the  particle  Stokes  number  increases  with 
gas  Mach  number  M|.  This  is  the  result  of  an  increase  in  the  relative  velocity  between  the 
particle  and  the  gas  with  Mach  number. 

Comparing  contours  of  T0  from  300  to  3500°K  for  7  =  1 .4  in  Fig.  8,  one  concludes  that 
an  increase  in  stagnation  temperature  (roughly  the  gas  temperature  back  of  the  shock) 
increases  the  gas  viscosity  and  thus  reduces  the  magnitude  of  the  particle  Stokes  number. 
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Figure  8.  Dimensionless  function  proportional  to  particle 
Stokes  number. 


4.4.2  Knudsen  Number 


The  particle  Knudsen  number  is  defined  as  the  ratio  of  the  mean  free  path  of  the  gas  to 
the  particle  diameter.  From  Eqs.  (10)  and  (11)  one  has 


Kn 


P2 

C2  G2  d 


(39) 


After  s.ome  rearrangement  using  the  procedure  of  the  previous  section,  one  obtains  an 
expression  in  terms  of  gas  stagnation  properties,  as  follows: 


where 


Kn  -  (t)“  (■*-)  f-’(M) 


(Tj/Ti) 

f*\ 

/  Po  \ 

[(T2/T,)  +  (TO/T,)(T0/TO)] 

l  P2  / 

V  Pl  / 

(40) 

(41) 


Equation  (40)  is  the  product  of  the  Knudsen  number  of  the  particle  in  the  stagnation  region 
and  a  dimensionless  function,  f2  (Mj). 
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Substituting  expressions  for  the  ratios  of  gas  properties  appearing  in  Eq.  (41)  from  Eqs. 
(1),  (2),  (28),  and  (29),  one  obtains  numerical  solutions  for  f2  (Mi),  and  these  results  are 
shown  in  Fig.  9.  As  indicated  in  Fig.  9,  f2  (Mi)  reaches  a  minimum  for  small  gas  Mach 
numbers  Mi  with 

For  constant  stagnation  conditions,  the  influence  of  gas  temperature  and  density  appears  to 
cause  a  decrease  in  the  Knudsen  number  behind  the  shock  at  small  Mi.  For  larger  Mach 
numbers  Mi,  the  gas  density  qj  decreases  rapidly,  accounting  for  the  increase  in  Kn. 


Figure  9.  Dimensionless  function  proportional  to  particle 
Knudsen  number. 
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4.4.3  Particle  Mach  Number 

The  particle  Mach  number  is  the  ratio  of  the  velocity  of  the  particle  relative  to  the 
ambient  fluid  to  the  local  speed  of  sound.  This  particle  parameter  is  summarized  in  Section 
3.2.1.  Equation  (5)  is  written  as 

MP2  =  f3  (Mj)  (42) 

and  numerical  solutions  for  this  expression  are  shown  in  Fig.  10  where  f3  (1)  =  0. 


Figure  10.  Particle  Mach  number  behind  a  normal 
shock  wave. 


4.4.4  Local  Reynolds  Number 


The  local  particle  Reynolds  number  represents  the  ratio  of  inertial  to  viscous  forces 
acting  on  the  particle  as  it  moves  relative  to  the  gas.  From  Eqs.  (7)  and  (11),  one  has  a 
maximum  value  of 


Re*  =  >  -  vjd  =  (~) 


1/2  Mp2 


Kn 
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where  Kn  and  MP1  are  defined  above  in  Sections  4.4.2  and  4.4.3  by  Eqs.  (40)  and  (42).  Thus, 
one  obtains,  in  terms  of  stagnation  conditions, 

Re2  =  Y1/2(— f4(Mj)  (43a) 

where 


f4<M,)  =  f3  (M,)/f2(M,) 


Numerical  solutions  to  Eq.  (43a)  in  Fig.  11  indicate  that  the  local  particle  Reynolds 
number  reaches  a  peak  at  M|  =*  2.  At  smaller  values  of  M]  ( <  2),  the  relative  velocity  of  the 
particle  with  respect  to  the  gas  behind  the  shock  wave  decreases  rapidly  and  f4  (])  =  0.  At 
larger  values  of  M]  ( >  2),  the  gas  viscosity  increases.  As  expected,  the  group  proportional  to 


Figure  11.  Dimensionless  function  proportional  to  local 
particle  Reynolds  number  behind  a  normal 
shock  wave. 
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the  particle  Knudsen  number  based  on  gas  properties  in  the  stagnation  region,  dgc/k, 
appears  in  Eq.  (43a)  as  it  does  in  Eq.  (40). 

4.4.5  Particle  Reynolds  Number 


A  particle  Reynolds  number  (Rep)  based  on  particle  velocity  ahead  of  the  shock  front 
and  on  gas  properties  behind  the  shock,  and  defined  by  Eq.  (22),  appears  in  the  normalized 
equation  of  particle  motion,  Eq.  (23).  The  magnitude  of  Rep  determines  whether  the  particle 
motion  can  be  approximated  as  Stokesian  or  non-Stokesian  and  thus  determines  the 
approximate  form  of  the  drag  coefficient  back  of  the  shock  front  as  indicated  in  Table  2. 


From  Eqs.  (22)  and  (43)  one  obtains 

Rep  =  Re2  ^  1  - 


or 


where 


(44) 

(44a) 


(44b) 


Numerical  solutions  to  Eq.  (44b)  are  shown  in  Fig.  12.  Similar  to  previous  calculations 
for  the  local  particle  Reynolds  number  shown  in  Fig.  11,  contours  of  f5  (MO  for  constant 
stagnation  temperature  peak  at  M|  «  2.  It  can  also  be  shown  that  f5  (1)  is  approximately 
equal  to  [(1  +  y)/2]l/l]  -t>. 


4.4.6  Similarity  Parameter 

When  the  particle  Reynolds  number  defined  by  Eqs.  (22)  and  (44a)  is  much  greater  than 
one  (i.e.,  Rep  >  102),  substitution  of  the  inertial  drag  coefficient,  Co,  from  Eq.  (19)  into  the 
normalized  expression  for  particle  motion,  Eq.  (23),  yields  the  similarity  parameter, 

St/Rep3.  Physically,  this  group  (like  the  Stokes  number)  represents  the  ratio  of  the  particle 
stopping  distance  to  the  probe  diameter.  Particle  trajectories  near  the  probe  inlet  will  be 
scaled  by  this  group  provided  Rep  is  large. 
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Figure  12.  Dimensionless  function  proportional  to  particle 
Reynolds  number,  Rep. 


Combining  Eqs.  (22)  and  (37),  one  obtains 


St 

Rep2/3 


(45) 


where 

Numerical  solutions  to  fe  (Mi)  show  an  increase  with  Mi  for  Mi  >  1.4  and  a  large 
dependence  on  the  ratio  of  specific  heats  as  shown  in  Fig.  13.  Moreover, 

f«0>  * 
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M1 

Figure  13.  Dimensionless  function  proportional  to  particle 
stopping  distance  for  non-Stokesian  particle 
(Re„  >  >  1) 

5.0  PARTICLE  BREAKUP 

Particle  agglomerates  and  liquid  droplets  are  subjected  to  large  shear  stresses  when  they 
encounter  shock  waves.  The  drag  forces  are  greatest  in  magnitude  at  the  instant  the  particle 
passes  behind  the  shock.  These  forces  decrease  in  magnitude  as  the  particle  equilibrates  with 
the  gas  downstream  from  the  shock.  The  problem  of  whether  agglomerates  and  droplets  can 
withstand  these  stresses  without  breaking  up  naturally  arises  in  many  types  of  particle  shock 
interactions. 

5.1  MAXIMUM  PARTICLE  DRAG 

The  drag  on  a  particle  is  given  by  Eq.  (13).  At  the  shock,  the  maximum  drag  is 

D  =  y  e2  Ap  CD  (V,  -  v2)2  (46) 
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We  now  consider  nozzle-generated  shocks  as  done  in  Section  4.  If  we  relate  gas 
properties  across  the  shock  and  into  the  stagnation  region,  Eq.  (46)  becomes 

D  -  T  Co  Ap  «.  Ci  MS  (-£-)  (fi)  (£)  (l  -  f)2  (47) 


Substituting  Eqs.  (1),  (28),  and  (29)  and  noting  that  e0  C2  =  P0,  one  obtains 


-4—  =  CD  g(M,) 
ApPo 


Since  g  —  0  for  M]  —  0  and  oo ,  we  set 

d  g(MQ  _  Q 
d  Mj 


(48) 


(49) 


(50) 


Solving  for  Mi  from  Eq.  (50),  one  obtains  for  the  position  of  the  maximum  in  g(M,), 

M,  =  (3  +  27) 1/2  (51) 

Substituting  Eq.  (51)  into  g(M|),  we  find  that  the  maximum  in  the  dimensionless  function  is 

*<">«)  -  -= - <52> 

[' +  (_LyJ“)(3  +  *0]  7“ 

Thus  if  7  =  1.4,  then  gmax  =  0.42  at  Mi  =  2.4.  Contours  of  g(Mj)  illustrating  these  features 
are  shown  in  Figs.  14  and  15. 

Numerical  solutions  are  now  sought  for  the  normalized  drag,  Eq.  (40),  with  constant  gas 
stagnation  conditions  and  particle  properties.  Since  the  drag  coefficient  is  of  the  form 


CD  =  CD  (Kn,Re2)  (53) 

where  the  arguments  appearing  in  Eq.  (53)  are  given  by  Eqs.  (40)  and  (43a),  values  of  the 
normalized  drag  were  determined  as  a  function  of  the  gas  Mach  number  Mj  and  the 
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parameter  de0/k.  These  results  are  shown  in  Fig.  14  for  air  (7  =  1.4)  and  in  Fig.  15  for  a  gas 
mixture  (7  =  1.16)  typical  of  solid  rocket  exhausts. 


Figure  14.  Dimensionless  function  proportional  to  particle 
drag  behind  normal  shock. 

The  results  presented  in  Figs.  15  and  16  demonstrate  that  particles  will  experience  the 
largest  stresses  at  a  gas  Mach  number  in  the  range  2  <  M!  <  2.5.  These  results  are 
insensitive  to  the  stagnation  temperature,  indicating  a  variation  of  less  than  5  percent  for 
stagnation  temperatures  in  the  range  500K  <  T0  <  3500K.  The  results  are  also  insensitive  to 
changes  in  Tq,  Sutherland’s  constant,  for  the  values  shown  in  Table  3  as  was  the  case  in 
Figs.  8  through  13. 

Expected  values  of  the  parameter  dg0/k  proportional  to  the  inverse  of  the  particle 
Knudsen  number  in  the  stagnation  region  are  indicated  in  Table  4. 
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Table  4.  Typical  Values  of  de0/k 


7  =  1.16 
PQ  “  600  psi 
d(/tm)  T0  =  3500°k 


7  =  1.4 
P0  =  50  psi 
T0  =  500°K 


1  50  -  28 


10  500 


280 


*R  =  2,86  x  106  CTn2-sec-2-‘>K- 1 
b  =  1.43  x  10~5  gm-cm- ’-see”  J-{0K) — 1/2 
k  =  .846  x  10“  ®  gm-cm-2 


5.2  STRENGTH  OF  AGGLOMERATES 


In  the  following  treatment,  we  assume  that  large  agglomerates  consist  of  small 
constituent  particles  that  are  uniform  spheres  of  size  dD.  An  example  of  an  agglomerate  is 
shown  in  Fig.  17.  Following  the  treatment  of  Rumpf  (Ref.  18),  we  also  assume  that  the  force 
of  adhesion  between  two  constituent  particles  is  F  and  that  each  small  sphere  of  size  dD  has  z 
contacts  (coordination  number).  If  the  voidage  is  1  -  <£,  the  number  of  spheres  per  unit 
cross-sectional  area  of  the  large  agglomerate  is 


n 


(54) 


Thus,  the  number  of  contacts  per  unit  cross-sectional  area  becomes 


nz  _  2<j>z 

2 


(55) 


d 


D 


Particle  Agglomerate 


Figure  17.  Schematic  of  particle  agglomerate. 
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where  the  factor  of  one-half  is  introduced  in  Eq.  (55)  since  each  contact  is  shared  by  two 
spheres. 


One  can  now  write  the  total  force  of  p  adhesion  per  unit  cross-sectional  area  of  the 
agglomerate  as 


P  _  0  z  F 
Ap  Tdfj 


(56) 


where  Ap  =  ird2/4  is  the  cross-sectional  area  of  the  agglomerate,  P  is  the  total  force  of 
adhesion,  and  F  is  the  particle  bonding  force  function.  In  Eq.  (56)  an  additional  factor  of 
1/2  has  been  introduced  by  Rumpf  (Ref.  18)  since  each  constituent  particle  must  break  only 
half  of  its  contacts  to  form  a  cleavage  plane.  Rumpf  also  demonstrates  that  z  (1  -  <t>)  »  3.1 
for  common  coordination  numbers  so  that  Eq.  (56)  simplifies  to 


P  _  d.  _F_ 
Ap  1  -  <b  d2, 


(57) 


It  now  remains  to  examine  the  various  forms  that  F  may  take. 


5.2.1  London-Van  der  Waals  Forces 


The  force  of  adhesion  between  clean,  dry,  uncharged  particles  and  surfaces  is  due  to 
London-Van  der  Waals  type  attractive  forces.  Nonpolar,  electrically  neutral  atoms  and 
molecules  have  momentary  dipoles  over  short  periods  of  time  due  to  the  movement  of 
orbital  electrons.  Integrating  the  attractive  forces  between  all  of  the  molecules  in  two  small 
spheres  separated  by  a  distance  “a”  gives  an  attractive  force  of  the  form 

F„  =  ^  (58) 

a*- 


The  distance  in  the  denominator  of  Eq.  (58)  is  restricted  to  a  <  2000A,  and  Bj  is  typically  in 
the  range  10-13  <  B|  <  10^ 11  ergs  (Ref.  14),  depending  on  the  nature  of  the  material  and 
the  media  separating  the  spheres. 


For  larger  separations,  a  >  2000A,  Lifschitz’  Eqs.  (18),  (14),  and  (19)  between  two 
spherical  particles  of  equal  diameter  are 


F 


v2  ~ 


B2do 


(59) 
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where  B2  *>  10-20  erg-cm  and  depends  on  particle  composition  and  the  media  separating  the 
spheres. 

It  should  be  noted  that  the  extensive  literature  available  on  the  subject  of  Van  der  Waals 
forces  demonstrates  that  the  calculations  suggested  above  are  only  approximate,  at  best,  as 
reviewed  by  Davies  (Ref.  14).  Electrostatic  charge  and  adsorbed  liquid  on  particle  surfaces 
can  result  in  considerable  error  in  the  simple  calculations  above. 

5.2.2  Electrostatic  Forces 


Suspended  particles  carry  electrostatic  charge.  There  are  many  mechanisms  of  charging, 
and  these  are  reviewed  by  Davies  (Ref.  14)  and  Friedlander  (Ref  8).  The  important 
mechanisms  fall  into  one  of  two  categories.  The  particles  are  charged  during  their  formation 
or  are  exposed  to  electrons  and  ions  carried  by  the  ambient  gas.  Assuming  two  spheres  with  a 
surface  charge  density  q,  separated  by  a  distance  "a”  much  smaller  than  the  diameter  of  the 
spheres  dD,  the  force  of  attraction  is 


Fe 


eir2  d£,q2 


(60) 


where  e  =  9  x  1018  dynes-cm2/couls  is  Coulomb’s  law  constant. 


In  the  case  of  particle  charging,  one  encounters  a  charge  limit  (Ref.  14).  This  represents 
the  maximum  electrical  charge  that  can  be  carried  by  the  particle  before  spontaneous 
emission  of  electrons  or  ions  occurs.  For  a  spherical  particle,  the  limiting  electron  or  ion 
charge  is  given  by  the  expression 


nP  ~ 


(61) 


where  np  represents  the  number  of  coulombs  of  charge  and  Kp  =  2,72  x  10-6  couls/cm2  or 
Kp  =  5.44  x  10-5  couls/cm2  for  the  electron  and  ion  limit,  respectively  (Ref.  14).  FromEq. 
(61),  and  noting  that 

np  =  rrd^q  (62) 


one  obtains  Kp  =  irq,^. 


Electrostatic  mechanisms  have  been  shown  experimentally  to  have  a  large  effect  on  the 
magnitude  of  the  force  of  adhesion  between  particles.  This  effect  is  much  larger  than 
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expected,  assuming  normal  charge  distributions,  suggesting  that  the  details  of  surface 
irregularities  and  non-uniform  charge  distributions  at  the  contacting  surface  are  not 
properly  accounted  for  in  existing  theories. 

5.2.3  Surface  Films 

When  particles  are  suspended  in  a  condensable  vapor,  a  film  of  liquid  develops  on  the 
surface.  If  the  vapor  pressure  is  near  its  saturated  value,  the  force  of  adhesion  between  two 
spheres  of  equal  diameter  with  smooth,  clean  surfaces  is  estimated  to  be 

Fs  =  2.2adD  (63) 

where  o  is  the  surface  tension  of  the  liquid  film.  In  practice,  the  force  of  adhesion  decreases 
sharply  below  the  value  given  by  Eq.  (63)  if  the  vapor  pressure  is  below  90  percent  of  its 
saturated  value.  Surface  contaminants  and  the  non-submergence  of  surface  irregularities 
will  also  reduce  the  expected  force  of  adhesion  (Ref.  14). 

5.2.4  Maximum  Strength 

The  previous  sections  have  outlined  a  number  of  important  forces  of  adhesion  which  are 
responsible  for  the  strength  of  an  agglomerate.  It  is  informative  to  consider  the  maximum 
force  for  each  of  the  mechanisms.  For  example,  in  the  case  of  Van  der  Waals  force  of 
attraction,  the  spheres  are  assumed  to  be  perfectly  smooth  with  a  minimum  distance  of 
separation.  This  would  correspond  to  roughly  10A  (equivalent  to  two  molecular  diameters) 
in  Eq.  (58)  and  2000  A  in  Eq.  (59).  The  remaining  parameters  necessary  are  the  surface 
tension  ( a  =  72  dynes/cm  for  water  at  standard  conditions)  and  the  electron  surface  charge 
density.  For  the  latter  value,  we  choose  20  electrons  as  the  average  charge  resulting  from 
static  electrification  on  a  1-^m  particle  (Ref.  19).  It  is  understood,  however,  that  the  total 
charge  could  increase  significantly  for  larger  particles  or  those  particles  exposed  to  large 
concentrations  of  unipolar  ions.  Substitution  of  these  values  into  Eqs.  (58),  (59),  (60),  and 
(63)  provides  expected  values  for  the  adhesion  force  on  a  1-nm  particle,  and  these  are  listed 
in  Table  5. 

It  is  apparent  from  Table  5  that  in  the  absence  of  a  limiting  charge  on  perfectly  smooth 
spheres,  either  Van  der  Waals,  electrostatic,  or  surface  tension  mechanisms  can  determine 
the  force  of  adhesion.  Exposure  of  the  spheres  to  large  concentrations  of  ions  of  a  single  sign 
(unipolar  charging)  may  produce  very  large  electrostatic  forces.  Moreover,  surface 
irregularities  may  alter  the  results  shown  in  Table  5.  Experimental  measurements  are 
necessary  to  resolve  the  question  in  specific  applications. 
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Table  5.  Maximum  Force  of  Adhesion  between  Two 
Smooth  Spheres  of  1-jtm  Diameter 


Mechanism 
Van  der  Waals 

Electrostatic 

Surface  film 


Force  (dynes) 

Fv,  =  10"2,  a  =  10  A 

Fv2  =  1.25  x  10“ 10,  a  =  2000  A 

Fe  =  2.4  x  10-2,  20  electrons 
a  =  10  A 

Fs  =  1.6  x  10-2,  7  =  72  dynes/cm 


5.3  STRENGTH  OF  ZnO  AGGLOMERATES 

Because  of  uncertainties  associated  with  the  simple  theories  of  particle  adhesion 
presented  in  Section  5.2,  experimental  data  are  usually  necessary  to  predict  the  conditions 
for  agglomerate  breakup.  In  this  section,  we  discuss  the  measured  strength  of  a  metal  oxide 
agglomerate  held  together  by  Van  der  Waals  forces.  These  results  may  be  applicable  to  the 
aluminum  oxide  agglomerate  known  to  exist  in  rocket  motor  plumes  (Ref.  3)  and  used  in  the 
seeding  of  the  gas  flow  in  high-speed  wind  tunnels  for  purposes  of  measurement  of  gas 
velocity  by  laser  doppler  velocimetry  (Ref.  20). 

Meissner  et  al.  (Ref.  20)  placed  ZnO  powder  of  constituent  particle  size  dD  =  0.13  and 
0.26  ftm  in  a  rotating  drum.  Mesh  sizes  in  the  range  from  0.25  to  1  mm  were  used,  and  the 
powder  was  tumbled  for  various  time  periods  to  produce  agglomerates  of  known  size  and 
volume  fraction  of  solids.  The  crushing  strength  of  the  powder  was  determined  by  placing 
the  material  between  microscope  slides  and  measuring  the  force  which  caused  failure,  as 
indicated  by  a  reduction  in  spacing  between  the  plates. 

The  results  of  Meissner’s  measurements  have  been  replotted  onto  a  single  curve  by  the 
appropriate  choice  of  the  function  represented  on  the  abscissa  of  Fig.  18.  Figure  18  indicates 
that  in  the  absence  of  large  amounts  of  water  vapor  (i.e.,  high  relative  humidity),  the 
strength  of  the  adhesive  force  F  within  the  ZnO  agglomerates  is  proportional  to  the 
constituent  particle  size  (dD).  This  suggests  that  Van  der  Waal’s  forces  are  the  mechanism  of 
adhesion  for  the  ZnO  agglomerates.  The  dependence  of  the  strength  of  the  agglomerate  on 
the  solid  fraction  (<£)  is  stronger,  however,  than  the  expected  form  presented  by  Rumpf  in 
Eq.  (57).  One  would  also  expect  the  force  of  failure  of  a  single  agglomerate  to  be  higher  than 
the  measured  value  in  Fig.  17,  since  a  void  existed  between  the  agglomerates  in  the  testing 
procedure  of  Meissner. 
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dQ 

Figure  18.  Crushing  strength  of  ZnO  agglomerates. 

5.4  STRENGTH  OF  LIQUID  DROPLETS 

Uniform  liquid  spherical  droplets  are  held  together  by  surface  tension  forces.  Since 
surface  tension  (rr)  has  units  of  force  per  unit  length,  the  strength  of  the  droplet  p  is  thus  the 
product  of  a  and  the  circumference  or 


p  =  cJTrd 


(64) 


Writing  Eq.  (64)  in  the  form  of  Eq.  (57),  one  obtains  a  strength  per  unit  cross-sectional  area 
of  the  droplet,  or 


(65) 


Caveny  and  Gany  (24)  have  studied  the  breakup  of  molten  AI/AI2O3  droplets  in 
accelerating  flow  fields.  Droplet  breakup  in  these  experiments  was  characterized  in  terms  of 
the  Weber  number  representing  the  ratio  of  inertial  to  surface  tension  forces.  Droplet 
distortion  is  typically  noted  at  a  Weber  number  of  4,  and  breakup  occurs  in  the  range  from 
12  to  20.  The  latter  values  were  slightly  smaller  than  the  values  measured  by  Ca%eny  and 
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Gany.  If  one  assumes  a  drag  coefficient  of  approximately  2.0  in  the  experiments  of  Caveny 
and  Gany,  for  the  particle  inertial  force  given  by  Eq.  (13),  droplet  distortion  occurs  when  the 
drag  per  unit  cross-sectional  area  of  the  droplet  equals  the  surface  tension  strength  given  by 
Eq.  (65).  Breakup  occurs,  however,  when  the  drag  force  is  about  four  times  the  droplet 
strength, 

5.5  CRITERIA  FOR  PARTICLE  AND  ZnO  BREAKUP 

From  the  foregoing,  we  are  now  in  a  position  to  establish  criteria  for  particle  breakup. 
The  discussion  is  restricted  to  the  possibility  of  breakup  behind  normal  shocks.  We 
hypothesize  that  the  onset  of  particle  breakup  occurs  when  the  drag  force  exceeds  either  the 
agglomerate  or  droplet  strength,  or  when  the  Weber  number  (Wb)  is  greater  than  unity. 

Wb  =  £  >  1  (66) 

where  D  is  given  by  Eq.  (49)  and  the  particle  strength  (P)  by  Eq.  (57)  or  (65).  In  Eq.  (57)  the 
constituent  particle  adhesive  force  (F)  which  appears  in  P  is  given  by  either  Eq.  (58),  (59), 
(60),  or  (63).  Equation  (66)  is  the  same  argument  used  to  predict  droplet  breakup  and  the 
definition  of  a  droplet  Weber  number. 

Substituting  for  particle  drag  and  strength  in  Eq.  (66),  the  following  criteria  are 
established  for  each  of  the  adhesive  mechanisms  discussed  in  Section  5.2.  The  criteria  for 
liquid  droplet  breakup  discussed  in  Section  5.4  is  also  included  for  completeness.  Thus, 
particle  breakup  is  expected  to  occur  for  the  conditions  listed  below. 

Van  der  Waals: 


Wb  =  dPp°cPg(Mi)  >  ] 
B  k,  (*) 


where 


Electrostatic: 


B  =  B,/a2,  a  <  2000 A 
B  =  B2/a3,  a  >  2 000 A 
k,(0)  =  0/(1  -  0) 


Wb 


P0CpG(M[)  >  t 
eq2k3  (0) 


(67) 


(68) 
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where 


kj  =  3.1xkj(<0) 


Surface  films: 


Wb 


dDP0CDg(Mi)  i 
■yMtf) 


(69) 


where 


k2(0)  =  2.17k,(0) 


Liquid  droplets: 


Wb 


dPpCpg(Mi) 

4cr 


(70) 


Equations  (67),  (68),  and  (69)  indicate  that  the  likelihood  of  particle  breakup  increases 
significantly  as  the  volume  fraction  of  solids  (0)  decreases.  Larger  gas  stagnation  pressures 
(P0)  or  constituent  particle  sizes  (dD)  also  contribute  to  agglomerate  breakup.  Moreover,  the 
product  CD  g(M|)  in  Figs.  15  and  16  increases  with  smaller  agglomerate  sizes  (d).  Thus,  we 
can  conclude  that  larger,  more  closely  packed  agglomerates  made  up  of  smaller  constituent 
particles  will  have  a  greater  tendency  to  resist  breakup  for  a  given  set  of  gas  stagnation 
conditions  and  shock  Mach  number  than  will  smaller  agglomerates  made  up  of  the  same  size 
particles. 

It  is  instructive  to  consider  the  results  of  the  experimental  measurements  of  the  strength 
of  ZnO  agglomerates  (Ref.  20).  Assuming  a  typical  range  of  values  of  2  x  10s  dynes/cm2  < 
h(d>)/dD  <  2  x  106  dynes/cm2  from  Fig.  18,  one  would  expect  breakup  to  occur  when  the 
drag  force  exceeds  these  values.  From  Table  4,  we  consider  an  agglomerate  of  size  d  =  10 

flowing  in  a  channel  with  a  stagnation  pressure  P0  =  50  psi  and  Mi  =  2.  Since  P0  =  3.4 
x  106  dynes/cm2,  and  dgQ/k  =  280,  andCpg(Mi)  =■  0.5  from  Fig.  15,  then  P0  CogCM])  = 
1.7  X  106.  Thus,  the  Weber  number  is  in  the  range  0.85  <  Wb  <  8.5,  and  breakup  may  be 
possible  in  this  case. 

If  the  same  10-^m  ZnO  agglomerate  were  suspended  in  a  gas  typical  of  rocket  plumes 
with  P0  =  600  psi,  7  =  1.16,  and  M]  =  4,  we  would  find  de0/k  =  5  X  102  from  Table  4  and 
calculate  P0  Cq  g(Mi)  =  8.1  x  106  dynes/cm2.  Thus,  the  Weber  number  is  in  the  range  2.0 
<  Wb  <  20,  and  breakup  is  much  more  likely  to  occur. 

It  should  be  noted,  however,  that  the  value  of  the  Weber  number  marking  the  onset  of 
particle  breakup  may  be  much  larger  than  unity,  as  in  the  experiments  of  Caveny  and  Gany 
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(Ref.  21).  Thus,  the  calculations  above  are  only  estimations  and,  of  course,  the  Weber 
number  would  be  considerably  reduced  in  the  transonic  or  hypersonic  regimes. 

6.0  SUMMARY  AND  RECOMMENDATIONS 

In  summary,  the  various  phenomena  which  influence  the  integrity  of  particle  samples 
obtained  from  blunt  probes  immersed  in  supersonic  flow  have  been  examined.  The  two  main 
sources  of  bias  were  identified  as  the  flow  around  the  probe,  which  is  complicated  because 
of  the  bow  shock  off  the  probe,  and  the  particle  breakup  possibly  induced  by  passage  of  the 
particle  agglomerates  or  droplets  through  the  bow  shock.  The  dimensionless  parameters 
necessary  for  scaling  these  phenomena  and  establishing  criteria  were  defined  (particle 
Reynolds,  Knudsen,  Mach,  Stokes,  and  Weber  numbers)  and  their  ranges  determined  for 
the  application  of  sampling  in  rocket  exhaust  plumes.  Expressions  for  particle  collection 
efficiency  were  defined  and  the  necessary  experimental  measurements  identified.  The 
breakup  of  agglomerated  particles  into  smaller  particles  was  examined  in  terms  of  the 
various  types  of  binding  forces.  The  only  material  for  which  such  binding  forces  have  been 
measured  is  zinc  oxide.  If  the  aluminum  oxide  from  rocket  exhausts  has  a  similar  binding 
force,  then  agglomerates  may  break  up. 

The  principal  benefit  of  this  study  is  the  definition  of  the  problems  to  be  solved  and 
identification  of  the  needed  data.  The  prediction  of  a  probe  collection  efficiency  is  much  too 
complicated  analytically  to  expect  believable  results;  experimental  measurements  are 
necessary.  The  experiments,  however,  need  not  be  over  the  entire  range  of  variables  since  the 
scaling  laws  have  been  established  by  this  study.  In  fact,  a  sub-scale  experiment  using  a 
laboratory  two-phase  flow  generator  and  probes  on  the  order  of  2-mm  diameter  rather  than 
the  25-mm  probes  used  in  rocket  exhausts  should  be  quite  acceptable.  Such  an  experiment 
can  be  used  to  quantify  most  probe  designs  for  collection  efficiency. 

The  second  area  requiring  experimentation  is  the  agglomerate  breakup.  It  is  first 
necessary  to  make  crushing  strength  experiments  on  aluminum  oxide  similar  to  those  made 
on  zinc  oxide.  In  this  experiment,  an  amount  of  material  collected  from  a  rocket  exhaust  is 
examined  for  size  distribution  and  then  placed  between  two  microscope  slides.  The 
measurement  consists  of  the  weight  required  to  produce  a  change  in  the  spacing.  This 
measurement  should  be  followed  by  an  experiment  in  which  the  effect  of  shock  waves  on 
particle  size  can  be  determined.  Either  a  free  jet  Mach  disc  or  a  shock  tube  would  be  used  for 
such  an  experiment. 

Finally,  it  should  be  noted  that  further  progress  in  this  area  is  hinged  on  experimental 
efforts. 
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NOMENCLATURE 

A  Cross-sectional  area  of  probe  (=  ttL2/4)  (cm2) 

Ae  Area  of  nozzle  exit  plane  (cm2) 

Ap  Cross-sectional  area  of  particle  (=  t d2/4)  (cm2) 

Ay  Area  of  nozzle  throat  (cm2) 
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a  Distance  of  separation  between  constituent  particles  (cm) 

B]  Van  der  Waals  constant  (ergs) 

B2  Van  der  Waals  constant  (erg-cm) 

b  Reference  value  in  Sutherland’s  equation  (gm-cm-l-sec~1-°K-1/2) 

C  Slip  correction  factor 

Cq  Drag  coefficient 

c  Speed  of  sound  (cm-sec-1) 

D  Particle  drag  behind  shock  (dynes) 

d  Particle  diameter  (cm) 

dD  Constituent  particle  diameter  (cm) 

F  Force  of  adhesion  (dynes) 

fi  Dimensionless  function  proportional  to  particle  Stokes  number 

Dimensionless  function  proportional  to  particle  Knudsen  number  behind  shock 
f3  Particle  Mach  number  ( =  MP2)  behind  shock 

fi  Dimensionless  function  proportional  to  local  particle  Reynolds  number 

behind  shock 

f5  Dimensionless  function  proportional  to  particle  Reynolds  number  Rep 

f6  Dimensionless  function  proportional  to  particle  stopping  distance  for  non- 

Stokesian  particles  (Rep  >  >  1) 

k  Constant  in  particle  groups  (  =  b/r1/2)  (gm-cm-2) 

kb  Boltzmann’s  constant  (gm-cm2-sec-2-°K_  ') 
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Kn  Particle  Knudsen  number  behind  shock 
Kp  Limit  of  surface  change  density  (couls-cm-2) 

L  Probe  diameter  (cm) 

M  Gas  Mach  number 

Mp  Particle  Mach  number 

m  Molecular  mass  (gm) 

mp  Particle  mass  (gm) 

n  Particles  per  unit  cross-sectional  area  (cm*2) 

ns  Particles  per  unit  mass  in  sample  probe  (gm-1) 

nm  Particles  per  unit  mass  in  mainstream  (gm-1) 

p  Particle  strength  (dynes) 

Pj  t2  Gas  pressures  (dynes-cm-2) 

q  Surface  change  density  (couls-cm-2) 

qm  Mass  flow  rate  in  mainstream  in  probe  area  (gm-sec-1) 

qs  Mass  flow  rate  in  probe  sample  line  (gm-sec  -  ‘) 

R  Specific  gas  constant  (cm2-sec-2-°K-  ■) 

Re  Local  particle  Reynolds  number  behind  shock 
Rep  Particle  Reynolds  number 

Rei,  Probe  Reynolds  number 

St  Particle  Stokes  number 
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T  Gas  temperature  (°K) 

Tg  Gas  Temperature  near  particle  (°K) 

Tp  Particle  temperature  (°K) 

To  Sutherland’s  constant  (°K) 

t  Time  (sec) 

u,  v,  w  Components  of  particle  and  fluid  velocities  in  three  orthogonal  directions 

(cm-sec  - ') 

v  Particle  velocity  (cm-sec  - !) 

v, -  Gas  velocity  (cm-sec  - 1) 

v+  Normalized  particle  velocity  ( =  v/V|) 

r*+  Normalized  gas  velocity  ( =  vf/vi) 

vr 

Wb  Particle  Weber  number 

z  Particle  coordination  number 

Greek  Symbols 

y  Ratio  of  specific  heats 

5  Shock  detachment  distance  (cm) 

e  Coulomb’s  law  constant  (dynes-cm2-coul- ') 

Q  Normalized  time  (  =  v'it/L) 

X  Gas  near  free  path  (cm) 

H  Gas  viscosity  (gm-cm- *-sec- ') 

q  Gas  density  (gm-cm-3) 
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gp  Particle  density  (gm-cm-3) 

a  Surface  tension  (dynes-cm- ') 

<j>  Solid  fraction  (agglomerate) 

a  Defined  function  of  Re 

Subscripts 

1  Upstream  of  shock 

2  Downstream  of  shock 

o  Stagnation  conditions 

r  Reference 

f  Fluid 
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